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We present exact results on the behavior of the thermodynamic Casimir force and the excess 
free energy in the framework of the d-dimensional spherical model with a power law long-range 
interaction decaying at large distances r as r~'*~'^, where a < d < 2a and < cr < 2. For a film 
geometry and under periodic boundary conditions we consider the behavior of these quantities near 
the bulk critical temperature T^, as well as for T > Tc and T < Tc. The universal finite-size scaling 
function governing the behavior of the force in the critical region is derived and its asymptotics 
are investigated. While in the critical and under critical region the force is of the order of L~'', 
for T > Tc it decays as L~'^~'^ , where L is the thickness of the film. We consider both the case of 
a finite system that has no phase transition of its own, when d — 1 < cr, as well as the case with 
d— 1 > cr, when one observes a dimensional crossover from d to a d— 1 dimensional critical behavior. 
The behavior of the force along the phase coexistence line for a magnetic field 11 = and T < Tc 
is also derived. We have proven analytically that the excess free energy is always negative and 
monotonically increasing function of T and H . For the Casimir force we have demonstrated that for 
any ct > 1 it is everywhere negative, i.e. an attraction between the surfaces bounding the system is 
to be observed. At T = Tc the force is an increasing function of T for cr > 1 and a decreasing one 
for CT < 1. For any d and ct the minimum of the force at T = Tc is always achieved at some H j^ 0. 
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I. INTRODUCTION 



When a fluid is confined in a film geometry with a 
thickness L, the boundary conditions which the order 
parameter has to fulfill at the surfaces bounding the sys- 
tem lead to a L dependence of the excess free energy. 
On its turn, the last lead to a force, conjugated to L, 
which is called the Casimir (solvation) force and the cor- 
responding effect - the thermodynamic Casimir effect. In 
this form it has been discussed for the first time by M. E. 
Fisher and de Gennes in 1978 [l|. The effect is dubbed 
so after the Dutch physicist Hendrik Casimir who first, 
in 1948 [21, predicted it considering the influence of the 
zero-point quantum mechanical vacuum fluctuations of 
the electromagnetic field on the resulting force between 
two infinite perfectly conducting planes placed against 
each other. In that form the effect is known as the quan- 
tum mechanical Casimir effect. For a long time the effect 
was considered as a theoretical curiosity but the inter- 
est in it has blossomed in the last decade. Numerous 
calculations and experiments have been performed both 
on the thermodynamic and the quantum Casimir effect. 
For a review on the thermodynamic effect the interested 
read_er_inight consult [j, |^ [^ , and for the quantum one 



The Casimir force in statistical-mechanical systems at 
a temperature T and in the presence of an external mag- 
netic field H is characterized by the excess free energy 
due to the finite-size contributions to the total free en- 
ergy of the system. In the case of a film geometry L x oo^ , 
and under given boundary conditions r imposed across 
the direction L, the Casimir force is defined as 



^Casimir i^, H , L) ~ - 



dff^{T,H,L) 
dL 



where f^^{T, H, L) is the excess free energy 
/r(T, H, L) = fAT, H, L) - Lfy^^i^T, H) 



(1.1) 



(1.2) 



Here /^ (T, H, L) is the full free energy per unit area and 
per fc^r, and fhuik{T, H) is the corresponding bulk free 
energy density. According to the standard finite-size scal- 
ing theory 0, ^^, under periodic boundary conditions 
T = p near the critical point T = Tc, H = (of the bulk 
system) one expects 

f^'^iT, H, L) = L-'^'^~^'ixf\atL^/\ hhL^/"), (1.3) 
wherefrom one has 



F, 



(p) 



Casimir 



{T,H,L) = L-'^41„,,(atLl/^6/.i^/'') 
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(1-4) 
Here the universal scaling functions of the free energy 

Xj'' {xi,X2) and the Casimir force -'^casimrX^i'^s) are 



related via the relation 

^climir(2;i,a;2) = 



{d-l)Xf\xi,X2) 

1 P\ 

Xl-z Xf{xi,X2) 

V oxi ■* 



A d 

'~^2^ 

V 0x2 



Ap) 



XY>{XUX2), (1.5) 



A and v are the standard critical exponents, a and h are 
nonuniversal metric factors, t = (T — Tc)/Tc is the re- 
duced temperature and h = /3H (/3 = (ksT)^^). We 
recall that, according to the general theory of the ther- 



(p) 

Casimir 



(xi,a;2) is 



modynamic Casimir effect y, |j, y], X, 
supposed to be negative under periodic boundary condi- 
tions (which corresponds to a mutual attraction of the 
'surfaces' bounding the system). The boundaries influ- 
ence the system to a depth given by the bulk correla- 
tion length ^ooiT) ^ \T — Td^" , where i/ is its critical 
exponent. When £,co(T) <^ L the Casimir force, as a 
fluctuation induced force between the plates, is negligi- 
ble. The force becomes long-ranged when £,oa(T) diverges 
near and below the bulk critical point Tc in an 0(n), 
n > 2 model system in the absence of an external mag- 
netic field [lllll2.ll3 |. Therefore in statistical-mechanical 
systems one can turn on and off the Casimir effect merely 
by changing, e.g., the temperature of the system. 

The temperature dependence of the Casimir force for 
two-dimensional systems has been investigated exactly 
only on the example of Ising strips |lj|. In 0{n) models 
for T > Tc the temperature dependence of the force has 
been considered in [ij . The only example where it is in- 
vestigated exactly as a function of both the temperature 
and of the magnetic field scaling variables is that of the 
three-dimensional spherical model with short ran ge i nter- 
action under periodic boundary conditions |l2l Il3l Il5l | . 
There results for the Casimir force in a mean-spherical 



model with L x 



geometry, 2 < d < 4, have been 



derived. The force is consistent with an attraction of the 
plates confining the system. In [ig some of the results of 
[13. [13 have been extended to a quantum version of the 
model. There the interaction has been taken to be long- 
ranged, with < cr < 2, where (7/2 < d < 3a/2, and the 
corresponding quantum phase transition has been con- 
sidered around T = 0. Very recently in [l5(, based on a 
derived there stress-tensor-like operator for critical lat- 
tice systems, the scaling functions of the force for the 
Sd Ising, XY and Heisenberg models have been obtained 
by Monte Carlo methods. The results suggest that, un- 
der periodic boundary conditions, the scaling function 

"'^Casimir (■^)/"' '^^ ^^^ ^^"^ 0{n) modcls practically coincide 
for large x, say, for x ~ L/S, ^ 2, where ^ is the true bulk 
correlation length. The last increases the helpfulness of 
the spherical model results (i.e. of the results in the limit 
n ^ 00), which are available in an explicit analytic form. 
Most of the results for the Casimir force are available 
only at T = Tc, i.e. for the Casimir amplitudes. They are 
obtained for d = 2 by using conformal-invariance meth- 
ods for a large class of models [3 . For d ^ 2 results for 



the amplitudes are available via field-theoretical renor- 
malization group theory in 4 — e dimensions [3 Hit Il7l | , 
Migdal-Kadanoff real-space rcnormalization group meth- 
ods [ia, and, by Monte Carlo methods [la, 113 . In addi- 
tion to the flat geometries some results about the Casimir 
amplitudes between spherical particles in a critical fluid 
have been derived too |17l l20l | . For the purposes of ex- 
perimental verification that type of geometry seems espe- 
cially suitable. For d = 3 the only exact ly k nown ampli- 
tude is that one for the spherical model [13. In the case 
d = a the amplitude is also known |l6l | for the quantum 
version of the model with long-ranged power-law interac- 
tion (in that case the amplitude in question characterizes 
the leading temperature correction to the ground state of 
the quantum system). 

It should be noted that in contrast to the quantum 
mechanical Casimir effect, that has been tested experi- 
mentally with high accuracy [23, [22, |23 |2J| (for a recent 
review on the existing experiments see, e.g. |25|), the 
statistical-mechanical Casimir effect lacks so far a quan- 
titatively satisfactory experimental verification. Never- 
theless, one has to stress that all the existing experiments 
|26l I27I l28l 1291 are in a qualitative agreement with the 
theoretical predictions. 

In this paper a theory of the scaling properties of the 
Casimir force of a spherical model with a power-law lead- 
ing long-ranged interactions (decreasing at long distances 
r as l/r'^+'^, with < cr < 2, and a < d < 2a) \s pre- 
sented. The results represent an extension to leading 
long-ranged interactions of the corresponding ones for 
system with short-ranged interaction |12. [13 ■ The lat- 
ter results, as we will see, can be reobtained by formally 
taking the limit ct — *■ 2~ in the expressions pertinent to 
the case of long-ranged interactions. 

All the interactions enter the exact expressions for the 
free energy only through their Fourier transform which 
leading asymptotic behavior is U{q) ~ a^rq" [1, l30l |. 
where a* — min(2, a). As it was shown for bulk systems 
by rcnormalization group arguments a > 2 corresponds 
to the case of sublcading long-ranged interactions, i.e. 
the universality class then does not depend on a p3l| and 
coincides with that one of systems with short-ranged in- 
teractions. Values satisfying < cr < 2 correspond to 
leading long-ranged interactions and the critical behav- 
ior depends then on a (see Ref. |33, 113 and references 
therein). In the current work we will restrict ourselves 
to the consideration of this case only. The other case 
of sublcading long-ranged interaction, i.e. when cr > 2 
is also of interest (involving, e.g., a serious modifica- 
tion of the standard finite-size scaling theory, see e.g. 
[33I |3J, I33, lla, I33), but will be considered elsewhere 

The investigation of the Casimir effect in a classical 
system with long-range interaction possesses some pecu- 
liarities in comparison with the short-range system. Due 
to the long-range character of the interaction there exists 
a natural attraction between the surfaces bounding the 
system. One easily can estimate that in systems with 



real boundaries (i.e. with no translational invariance) in 
the ordered state the L-dependent part of the excess free 
energy that is raised by the direct inter-particle (spin) in- 
teraction is of the order of L~'^~^^. In the critical region 
one still has some effects stemming from that interac- 
tion on the background of which develops the fluctuating 
induced new attraction between the surfaces that is in 
fact the critical Casimir force. In the definition (|1.1|) 
used here, that is the common one when one considers 
short-range systems, these effects are superposed simul- 
taneously. In the current article we will investigate their 
interplay. An interesting case when forces of similar ori- 
gin are acting simultaneously is that one of the wetting 
when the wetting layer is nearly critical and intrudes be- 
tween two noncritical phases if one takes into account 
the effect of long-range correlations and that one of the 
long-range Van dcr Waals forces |39l liol l4ll| . 

The structure of the article is as follows. In Sec- 
tion ni wc briefly describe the spherical model (which, 
in systems with a translational invariance, is equivalent 
to the n ^ oo limit of the 0{n) models) and give all ba- 
sic expressions needed to investigate the behavior of the 
Casimir force. In Section IlIII we derive the scaling func- 
tion of the excess free energy and the Casimir force, and 
investigate the leading asymptotic behavior of the force 
both above and below the critical point. In Section lTvl we 
consider in some details the behavior of the force along 
the phase coexistence line T < Tc, H ~ 0. In Section Ivl 
we investigate the monotonicity properties of the excess 
free energy, and the Casimir force, and prove analytically 
that both the excess free energy and the force are nega- 
tive for any T and H (for cr > 1). The last implies that 
the force between the boundary surfaces of the system 
is always attractive. The article closes with a discussion 
given in Section IVTI The technical details needed in the 
main text are organized in a series of Appendices. 



II. THE MODEL 

We consider the ferromagnetic mean spherical model 
with long-range interaction confined to a fully finite d- 
dimcnsional hypercubic lattice Cd of N = \Cd\ sites. 
The model is defined by 



n 



- 2_^ Jij StSj - Hy tSi , 



(2.1) 



where Si is the spin variable at site i, Jij is the interac- 
tion matrix between spins at sites i and j, and H is an 
ordering external magnetic field. The long-wave length 
asymptotic form of the Fourier transform J'(q) of the 
interaction potential Jij is 

j(s^^j{q)\\-p,ujM)\. |q|^o, p, >o. 



for large spin separations r ~ |r|. The spins in the model 
under consideration obey the spherical constraint 



Y.^si 



N, 



(2.2) 



where (• • • ) denotes standard thermodynamic averages 
taken with the Hamiltonian Ti. and N is the total number 
of spins located at sites i of finite hypercubic lattice Cd 
of size Li X L2 X ■ ■ ■ Ld = N (here Li are the linear sizes 
of the system measured in units of the lattice constants). 
Under periodic boundary conditions imposed along the 
finite directions of the system, the free energy density of 
the model is given by 



fi 



/3Td,A0.H,-L\A) = -sup{t/rf,,(0,L|A) 



0>O 



PJ{0)Po 

2tt 



PH^ 



J{'i)Pa(^ 



In 



-/3J(0)p, (<^+— ) }, (2.3a) 



where 



[/,,,(^,L|A) = l^ln[^ + c.<,(q)]. 



(2.3b) 



Here the vector q has the components {51,(72,' ■■ , 9d} 
where qj = 2Tmj / Lj and Uj G {—Mj, • • • , Mj — 1} with 
Mj = LjAj/{2tt) :^ 1 being integer numbers, and A^ the 
cutoff at the boundaries of the first Brillouin zone along 
the J direction. The spherical field is introduced to 
ensure the fulfillment of the constraint H2.2|l . It is the 
solution of the equation 



PJiO)pa 1 - 



H^ 






<fJ^{Q)pl) N^cf) + uM) 



(2.3c) 
Equations (|2.3a|) and H2.3c|) contain all the necessary 
information for the investigation of the critical behavior 
of the model under consideration. 

In the bulk limit, when all the sizes of the system are 
infinite, the d-dimcnsional sums over the momentum vec- 
tor q in Eqs. Ij2.3bp and (|2.3c|l transform into integrals 
over the first Brillouin zone. For example one has 



UdA'P\^) 



1 



dqi 



\Ti[(l) + uj„{qi,q2,- 



dqd 



-A 



,qd)]. (2.4) 



By analyzing the equation for the spherical field H2.3c|) 
in the bulk limit it is easy to show that the system ex- 
hibits a phase transition for d > a at the critical point, 
Pc, given by 



We suppose that the interaction in the system is long- 
rangcd with < cr < 2, i.e. aJo-(q) — Iq^. This corre- 
sponds to the inverse power-law behavior J'(r) ~ r~'^~'^ , 



f3cJ{0)pa 



(27r)« 



dqi 



dqd — -, r- 

A ^<j(qi,q2,--- ,qd) 
(2.5) 



III. SCALING FORM OF THE EXCESS FREE 

ENERGY AND THE CRITICAL CASIMIR 

FORCE 

In the remainder we consider a system with a film ge- 
ometry L X 00*^^^, which resuhs after taking the hmits 
L2 ^ oo, • • • , L^i — > 00 and setting Li = L. For the sim- 
plicity of notations we will only consider the case when 
all cut-off variables are taken to be equal to each other, 
i.e. A; = A, z = 1, • • • , d. Then Ud,tT{<l>i L|A) becomes 



f/d.<.(</),L|A) 



dq: 



'2 • • • 



dqd 



ln[(f) + uj^{qi,q2, • • • , qd)]- (3.1) 



The above sum can be evaluated using the Poisson sum- 
mation formula and the identity 

r°° dr 
ln(l + z")=a/ — (l-e-^^O^a(-x''), (3.2) 
Jo 2; 

where Ea{x) = Ea^i{x)^ and 



(3.3) 



are the Mittag-Lefflcr functions. For a review on the 
properties of Eafj{z) and other related to them func- 
tions, as well as for their application in statistical and 
continuum mechanics see Ref. [i^ (see also Ref. |4J|). 
The properties used in the current article are summa- 
rized in Appendix 1X1 

After some algebra for the full free energy density we 
receive: 

l3TdAP. H, L) = pj'dAP. H) - ^L-'^ICdAL'^^), (3.4) 



where 



and 



J^dAP,H)= lim TdAP,H,L), 

L—^oo 



00 pQQ 



/2-1 



exp 



4a; 



x£;./2,i i-x'^/^y 



(3.5) 



The main advantage of the above expression for the free 
energy, despite its complicated form in comparison to 
equation H2.3a(l . is the simplified dependence on the size 
L which now enters only via the arguments of some func- 
tions. This gives the possibility, as it is explained below, 
to obtain the scaling functions of the excess free energy 
and the Casimir force. It is worthwhile noting that under 
a sharp cutoff A a special care has to be taken when per- 
forming finite-size scaling calculations in order to avoid 



receiving artificial, i.e. not existing in real systems, finite- 
size A-dependent contributions. This question is consid- 
ered in details in |23|. In obtaining H3.4() the suggested 
there receipt has been applied (see Eq. (27) in [33 and 
the discussion devoted to it). According to these find- 
ings, for the finite-size contributions in the following we 
are going to send the cutoff to infinity. 

In equation H3.4|) , cj) is the solution of the correspond- 
ing spherical field equation that follows by requiring the 
partial derivative of the right hand side of equation l|3.4l) 
with respect to (j> to be zero. Let us denote the solution of 
the corresponding bulk spherical equation by (f)oo ■ Then, 
for the excess free energy (per unit area) it is possible 
to obtain from equations H1.2() and H3.4|l the finite size 
scaling form, valid for a < d < 2ct, 



n/], H, L\d) = l3L~^^-^^Xf (xi, X2) 
with scaling variables 

xi = (/3-/3,)J(0)p,Li/^ 
and 



X2 = iJL^/ V/3/^(0)Pa 



(3.6) 



(3.7a) 



(3.7b) 



Here v = \/{d - a) and A = (d -h a)/{2{d - a)) are the 
critical exponents of the spherical model (for a < d < 2a, 
and < cr < 2). In equation (|3.6|l the universal scaling 
function X°^ (xi,a;2) of the excess free energy has the 
form 



Xf{xi,X2) 


1 2^ M 1 

2 ^\yL yooj 2 


i{yL- 


-yoo) 


-^\DdA{yl^^~yi/>'')-li^<iAyL), 


(3.8) 


where the y^ ^ 4>lL" , yoo ^ (t'ooL" ■ and 




Dd.a = 2^ 


(ATrf/^T - cr sin — 


-1 


(3.9) 



In Eq. H3.8|l j/l is the solution of the spherical field 
equation for the finite system obtained by minimizing 
the free energy with respect to yL 



yi 



d/a-l d_ 

dyL 



iCdAyL). (3.10) 



For the infinite system the corresponding equation is 



a;i = -5 — \DdAy'. 



d/cr-l 



(3.11) 



According to equation (|1.1() . the finite-size scaling 
function of the Casimir force for the system under con- 
sideration is 



-'fCasimir (3^1,3^2) 



cr+1 , / 1 1 



yL yoo 



Cr-1 



xi {yL - yoo) 



-(^-1),. .f,f^_,i/^) 



2d 



'-\DdA {l 



-{d~ l)lCd,a{yL)- 



(3.12) 



Note that in the hmit cr ^ 2" Eqs. (|^ - (|X1H) 
reproduce exactly the correspondin g on es for the case 
of the short-range interaction [l2, ll^l ll^ . In such 
a case the above equations simphfy greatly since then 
Ei^i{z) = exp(z), and the function ICd,a{y) defined in 
Eq. (|3.5|l becomes 



l^dAy) 



(27r)''/: 



;y 



d/4 



Y.^-"^Ka/^{l^), (3.13) 



1=1 



where K^ is the modified Bessel function. 

In the present article wc will concentrate on the in- 
vestigation of the behavior of the Casiniir force and the 
excess free energy in different regions of the phase dia- 
gram. We will also evaluate some critical amplitudes for 
selected values of the parameters d and a. The analysis 
will be done analytically for the cases where one can ob- 
tain simple expressions and is then extended numerically 
to cases which are not accessible by analytical means. 

In FIG. n we present the numerical evaluation of the 
Casimir amplitudes as a function of d for some selected 
values of a. The results show that the amplitude is an 
increasing function of d at fixed a, and an increasing func- 
tion of (T at a fixed d. Note also that in accordance with 
the general expectations, the amplitudes are negative. 
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FIG. 1: Behavior of the Casimir amplitude as a function of d. 

In order to obtain the amplitudes, one needs to know 
the value of yiiT) at the critical point T = Tc that is 
the solution of the equation for the spherical field H3.1U|I . 
These results have their own important physical meaning. 
We recall that yi, is directly connected to the finite-size 
correlation length ^^ = Lyj^ of the system Q- The 
results for yL{Tc) are shown in FIG.|21 

In FIG. 121 we present our results for the Casimir force 
evaluated at the bulk critical point of the model as a func- 
tion of d for some selected values of a. We observe that 
the Casimir force behaves in a different way depending 
on whether <t is smaller or larger than a = 1. For cr < 1 
it is decreasing monotonically as a function of d, while 
for (7 > 1 it is not. 
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FIG. 2: Behavior of the scaling variable ijl as a function of 
d at the critical point T = Tc. We recall that the finite-size 



correlation length ^l is related to yL via ^l = Ly] 
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FIG. 3: The behavior of the Casimir force at T 
function of d. 
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In the following wc turn our attention to the inves- 
tigation of the thermodynamic functions of interest as 
a function of the scaling variable xi for fixed d and a. 
Let us first consider the situations where it is possible to 
obtain sonre results analytically. 

Let us first consider the asymptotic of the excess free 
energy and the Casimir force in the under critical region 
(i.e. T < Tc). Taking into account that then (i.e. when 
xi > 1, a;2 = 0), according to Eqs. H3.in() an d (|3.11() 
i/L -^ 0+, yoo ~ 0, as well as the asymptotic (jB13|l of 
the function ICd.aiVL) for snrall values of the argunrent 
(derived in Appendix Ib)| it is easy to see that below the 
critical temperature 



Xf (xi -^ c«, 0) ~ — 



27r'^/2 



i) cid), 



(3.14) 



and 



X 



Casimir {X, -. OO, 0) ^ "^^^^ (^ ) ^('^)- ^^-IS) 



The above results reflect the dominating contribution of 
the Goldstone modes in the under critical-regime of an 
0{n) model - both the excess free energy and the Casimir 
force do not tend exponentially-in-i to zero, but to finite 
constants. For a = 2 these constants coincide with those 
known from short-range systems (see, e.g. [13 and the 
references cited therein). Note also that, in contrast with 
systems with real boundaries, the direct inter-spin long 
ranged interaction below Tc does not lead to a L~('^~^' 
contribution, which is well known from studies of van der 
Waals systems exhibiting wetting phase transitions [33, 
|4J| . This is due to the application of periodic boundary 
conditions, i.e. the system under consideration lacks real 
physical boundaries. 

Let us consider the critical behavior of the force for T > 
Tc in a bit more details. Then, when a;2 = and xi -^ 
— cx) from (|3.1U|) and H3.11|) one obtains yi ~ yoo{^+£d.a), 
where 



£d,, 



ad,i 



I n I d/a+l ' 



and 



2/0 



\xi\ 



\D, 



(3.16) 



(3.17) 



Therefore, the leading behavior of the scaling function of 
the force in that region is 



X, 



Casimir 



where 



-Ad,ayJ 

-Ad,a [iPc - P)J{0)p^/\Dd, 



A,,. = ^(a + d-l). 



(3.18) 



(3.19) 



Eq. (|3.18|) implies that above T^ -Fcasimir — 
-X+\t\-^L-'^'^+''\ with 7 = (T/(d-cr), andX+ > 0, i.e. 
the force remains attractive and decays in a power-in-i 
and not in an exponentially-in-L way, as it is in systems 
with short ranged interactions. This behavior is in full 
correspondence with the long-ranged character of the in- 
teraction. Similar is, as it has been recently established, 
also the behavior of the Casimir force and the excess free 
energy in systems with van der Waals type interaction 
[37I I (see also |23), despite that their critical exponents 
are that ones of the short-ranged systems. 

The obtained analytical results are supported by nu- 
merical analysis of the expressions for the scaling func- 
tions of the excess free energy and the Casimir force at 
zero external field. The corresponding data is presented 
in FIG. 01 (for the excess free energy) and in FIG. [S] (for 



the Casimir force). While the scaling function of the 
excess free energy is monotonic regardless of the used 
values of d and ct, the behavior of Casimir force depends 
strongly on the range of the interaction a. For ct > 1 
it is monotonically increasing as it can be seen from the 
case (7 = 2, corresponding to short range interaction, and 
the long-range case with a = 1.5. For ct = 1 the mono- 
tonicity changes and Xcasimii(a;i,0) becomes decreasing 
for values of ct < 1. As example we show its behavior for 
CT = 0.75. 
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FIG. 4: The universal finite-size scaling function of the excess 
free energy Xf{xi,Q) from Eq. 13.811 as a function of —x\ ~ 
(T — Tc)L^'" , for some selected values of a at zero external 
magnetic field. One observes that, in full accordance with 
the corresponding statement from Section [3 ^Casimir (a;i,0) 
is a monotonically increasing function of the temperature T 
regardless of value of o. 

We close this section by presenting the outcome of the 
numerical analysis of the behavior of the scaling func- 
tions of the excess free energy, shown in FIG. [HI and that 
of the Casimir force, shown in FIG. |7| as a function of 
the scaling variable X2 at the bulk critical temperature. 
One observes that the excess free energy is a monoton- 
ically increasing function of the external magnetic field 
H independently of the range of the interaction. How- 
ever the Casimr force is a nonmonotonic function of H 
and has a minimum at X2 ^ which depth depends of 
the parameter ct. The minimum is found to be at X2 ~ 
0.084, 0.145, 0.263 and 0.416 for ct = 2, 1.5, 1 and 0.75, 
respectively. So, as long as ct goes smaller the minimum 
becomes deeper. Indeed the ratio of the Casimir force 
evaluated at the minimum to its value at H = is a de- 
creasing function of ct. It is given by 1.017, 1.073, 1.215 
and 1.513 for ct = 2, 1.5, 1 and 0.75, respectively. 



IV. CASIMIR (SOLVATION) FORCE ALONG 
THE PHASE COEXISTENCE LINE 



Here we investigate the behavior of the Casimir force 
along the line H = when T < Tc- This is a line of a first 
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FIG. 5: The universal finite-size scaling function of the 
Casimir force Xcasimir(a;i,0) as a function of the scaling vari- 
able —xi ~ (T — Tc)!/^'", at zero external magnetic field 
H = 0. One observes that, in full accordance with the corre- 
sponding statement from Section[Vl Xcasimir(2;i, 0) is a mono- 
tonically increasing function of the temperature T (for ct > 1) 
and possesses a complex behavior for ct < 1. 



o" 



O 



0.0 



-0.1 



-0.2 



-0.3 





' :::"::::::':-"■■-■ -.-.-".-J^ 


/f -0.1 

:i -0.2 

fl -0.3 


^^.i,,..--— — 




1 




2 4 


; 0=1.00, « 
; a =1.50, a 


= 1.50 

= 2.00 


1 


1 


1 



20 40 60 

Scaling variable x. 



80 



FIG. 7: The universal finite-size scaling function of the 
Casimir force Xcasimir, for some values of a, as a function 
of the scaling variable X2 ~ HL '^ at the bulk critical tem- 
perature r = Tc- One observes that, Xcasimir (0,0:2) is not a 
monotonically increasing function of the field H for all values 
of a < 2 including the short-range case. 
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FIG. 6: The universal finite-size scaling function of the excess 
free energy Xf (0,2:2), for some values of a, as a function of 
the scaling variable X2 ~ HL '" at the bulk critical point 
T = Tc- One observes Xf{0, X2) is a monotonically increasing 
function of the field H for arbitrary a. 



order phase transition with respect to the magnetic field 
H. The finite-size rounding of the first-order transitions 
in 0{n) models has been already studied by Fisher and 
Privman in [43 for a fully finite and cylinder geometries. 
Later their predictions have been verified in details for 
the spherical model system with such a geometry in jig , 
while in |47l liq their arguments have been extended to 
a geometry of the type L'^~'^ x oo'^ , where d' has been 
chosen so that no phase transition of its own exists in the 
finite system, i.e. d' < a has been supposed. Here we 
extend these investigations to cover also the cases d' = a 



and d' > a in systems with a film geometry, i.e. when 
d' = d — 1. We will be only interested in the behavior of 
the Casimir force. 

For T < Tc and small H Eqs. (|T^ - (|TT^ are still 
valid, but there the limit j/_l <C 1 has to be taken (i.e. we 
suppose that xi ^ x"^)- As it is clear from Eq. IJB13|I . 
then there are three subcases to be considered. 

i) The case d — 1 < a. 

Then in the finite system there is no phase transition 
on its own. For the excess free energy one obtains 
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(4.1) 



(4.2) 



Here niL = H/[pcrJ^{Q)(t)] is the magnetization of the fi- 
nite system, m,Q = ^\ — T jTc is the spontaneous mag- 
netization, and x^n ~ (3mo{T)L£^j^~ H, which has the 
meaning of the ratio of the total magnetic energy in the 
correlated volume 14or = L£^l~^ to the thermal energy 
ksT per degree of freedom, is the scaling variable. (We 
recall that in the spherical model the true finite-size cor- 
relation length ^L is equal to 4''^/'^ IE E3-) Next, it 
is easy to see from Eq. (|4.2|l that Xm = 0(1) involves 
H = 0(L-'^/(i+'^-'^)), that is the scale on which the 
jump in the bulk magnetization is rounded off. From 



this observation and from Eq. H4.1|l one obtains that the 
H dependent correction to the Casimir force is then of 
the order of L^'^/i^+'^-d) ^j^^^^ ^Yi&t cr/(l + a - d) > d 
for d > a, and, so, the term proportional to H in Eq. 
(|4.1() will indeed contribute as a correction towards the 
Casimir force). 

a) The case d — 1 = a. 

This is the borderline case between that one when in 
the finite system there is no phase transition of its own 
(for d—1 < a) and that one in which in the finite system 
there is such a phase transition (for d — 1 > a). In this 
case an essential singular point exists in the finite-size 
system a.t T = H = 0. For the excess free energy one 
now obtains 



rif3,H) 



2T^dl2 



T{d/2)Q{d)L 



-(d-i) 



+(3moHL ll-^ 
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(4.3) 
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(4^)-/2r(a/2) x™ 
(4.4) 
and .T,„ = l3mo{T)HL^'[-^ / \n{L/^L)- The above equa- 
tions are to be compared with the previous case. One 
observes, that the main difference is the existence of 
logarithmic-in-L dependence that is introduced via the 
scaling field variable a;,„. As a result the rounding of the 
jump in the magnetization takes place on a scale given 
hy H = L~'^ exp(— const. L), i.e. the scale in this case is 
exponentially small in L. 

ii) The case d ~ 1 > a. 

In this case there is a true phase transition of its own 
in the finite system at some Tc^l ~ Tc — eL^^^'^, i.e. no 
rounding of the jump of the magnetization is possible. 
One only observes i-dependent corrections of the finite- 
size magnetization rriL with respect to the spontaneous 
magnetization toq. One finds that the crossover from d 
to d — 1 critical behavior happens at Tc^l with 



^ (27r)- T{d/2) l3cJ{Q)p.' 
and, when \H\L" < 1, 



(4.5) 



rif3,H) = --^Tid/2)C{d)L-^''-'^+f3moHLj^, 

(4.6) 

(4.7) 



with 



r^'~'^/' C,.a 



2(27r)- T{dl2)(3ml{T)J{Q)p,' 



and rriL — to,o(1 — a/2). 

Finally, we would like to note that in 0{n) systems 
one observes for T < Tc in addition to the rounding of 
the jump of the order parameter also rounding of the 



spin wave singularities. According to the general the- 
ory [4J, |4g , their rounding occurs on the scale for which 
Xs = \H\L'' ^ 0{1). As it is clear from Eq. ^J^ (and 
taking into account that if T < Tc one can rewrite xi as 
xi = f3mo{T)^pcJ^{0), with xi ^ 1) the scale on which 
the rounding of the spin wave singularities sets in involves 
that xi ~ 2^2 there. Then, in this regime, the solution 
of the spherical field equations for the finite and the in- 
finite system (I3.10|) and (|3.11() will be again y^ = 0{1) 
and yoo — 0(1). Since xi and X2 can be expressed from 
Eqs. H3.10|) and (|3.11|) in terms of yL and j/tx., we con- 
clude that, according to (|3.12|) . in the regime in which 
the spin waves are of importance, the Casimir force will 
be Fcasimir = 0(L^'^), posscssing a nontrivial H depen- 
dence. If one would like to reveal more on this depen- 
dence the numerical treatment is unavoidable. Note that 
when the field is strong enough to suppress the spin- wave 
excitations, i.e. when a;^ 3> 1 and T < T^ one will have 
an Ising-like system. In this regime y^ ^ 1, j/oo ^ 1, 
and the Casimir force will be of the order of L^^^+cr) j-ggg 
Eq. (|3.18|) ') under periodic boundary conditions. (If the 
system was possessing real bounding surfaces like, say, 
under Dirichlet-Drichlet boundary conditions, one would 
expect that the corresponding contribution in the force 
is of the order of L^°'.) 



V. MONOTONICITY PROPERTIES OF THE 

EXCESS FREE ENERGY AND THE CASIMIR 

FORCE 

Let us denote by g^ {x2 , y) and goo {x2 , y) the right- 
hand side of Eqs. H3.10|l and H3.11|l . respectively. Now 
we prove that 

i) gL{x2,y) > 9oo{x2,y) and 

ii) that gL(x2,y) and gooix2,y) are monotonically de- 
creasing functions of y. 

i) First, let us note that Ea^—x) is a completely 
monotonic function of a; > p9ll50ll5lll53 | for < a < 1 
and (3 > a. (In [43 this property was shown to hold for 
EnA j—x) = En{—x) and was later extended to Ea^i3{—x) 
in |50| and [51|; see also [521 •) This means that for all 
n = 0, 1, 2, 3, • • • one has 

(„^)n ^"£^a/3(-x) > 0, :^ > 0, < a < 1, /? > a. (5.1) 
ax" 

Then, from n = it immediately follows that 
Ea,a{—x) > when a; > 0. Now, from Eqs. H3.10|l and 
(I3.11|) . it immediately follows that (7^(2:2, y) > 5oo(a^2, J/)- 

ii) The required property follows from the monotonic- 
ity of the function Ea^ai—x) for x > and the explicit 
form of the right hand sides of Eqs. H3.1U|I and (|3.11|) . 

Having proved i) and ii), it is easy to understand now 
that for any given values xi and X2 of the scaling variables 
the solution of the spherical field equation for the finite 
system will be larger than that for the infinite system, i.e. 
j/l(xi,X2) > yoo{xi,X2)- (Since the correlation lengths 
in the finite and the infinite system are ^l = Vl ^-^^ 
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FIG. 8: The universal finite-size scaling function of the Casimir force as a function of scaling variables xi and X2 for some 
values of the parameter a and the corresponding values of d. The visualization is limited to positive values of X2 since the 
function is even in H. 



Coo = 2/oo 13 J correspondingly, the physical meaning 
of the above result is that the correlation length of the 
finite system is always smaller than that of the infinite 
one.) We are then ready to prove that 

A) For Xi > and 2:2 = the excess free energy scaling 
function is negative, i.e. Xf(xi > 0,a;2 = 0) < 0. 

B) The excess free energy scaling function Xf{xi,X2) 
is a monotonically increasing function of the temperature 
T and the magnetic field \H\. 

Let us start with statement A). 

A) From the explicit form of the Eq. H3.8|l it is clear 
that the statement A) will be true if Ea,i{—x) > when 
X > 0. The last follows from (|5.1|) for n = 0, and, thus, 

Xf{xi,X2) <0. 

Let us now prove the statement B. 

B) Prom Eq. H3.8|) one obtains 



dX 



f 



1 



and 



dxi 



dx9. 



(2/00 - Vl) < 0, 



^ 2;2( — 
2/00 



VL 



-)>0. 



(5.2) 



(5.3) 



Eq. H5.2fl implies that Xf{xi,X2) is a monotonically in- 
creasing function of T, whereas Eq. H5.3|l states that it 
is a monotonically increasing function of \H\ too. 



Using B) one can now prove that: 

C) The excess free energy scaling function is negative 
for any T and H , i.e. Xf{xi, X2) < for any xi and X2. 

Indeed, from the monotonicity property B) and from 
A) it is clear that in order to prove C) it is enough to 
show that it holds for values of T above Tc, i.e. when 
j/L > 1 and j/oo > 1. T hen, from Eqs. ifXTUIl . (|TTT|) 
and the asymptotic (|B14|) one obtains y^ = yoo(l + £), 
< e <C 1, where 



a-d,, 



\D,,,\y'ol^~\^-i) 



2^ 



(5.4) 



Next, from Eq. (|?^ it follows that 

Xf{x^,X2):^-l^^{l-e)<0. (5.5) 

^ 2/00 

Thus, the excess free energy is indeed always negative. 

Finally, we prove that 

D) For a > 1 the Casimir force is always negative, i.e. 
it is a force of attraction between the surfaces bounding 
the system. 

We start by multiplying Eq. H3.11|l with y^o and Eq. 
(I3.1U|) with y^, and then adding the results together. One 
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obtains 



11^ ^ \ in, I (.Al'^ _,,dl<y 



UL J/oo 

VL-, ICd.aiVL)- 

dVL 



(5.6) 



Inserting the above expression in Eq. H3.12|l . one obtains 



-^Casimir (a;i,a;2) = xl 



1 1 

VL 2/oc 



-\ (i - ^) \D,A {y'l!" - y'^^) - \{d - i)K.aAyL) 



2 dyL 



(5.7) 



Since, according to what aheady has been proven, y^ > 
yoc, and JCd.aiyh) is a positive and monotonically de- 
creasing function of y^ (the last follows from the explicit 
form of 1Cd,a{yL) given in Eq. (|3.5|) and the property 
(|5.1() of Ea^x) for n = and n = 1), from the above 
expression one immediately confirms the validity of state- 
ment D). In addition, from Eq. (|3.12(l it is easy to see 
that Xcasimir (a;i,a;2) < also for cr < 1 if xi < 0, i.e. 
for T > Tc. Furthermore, from Eqs. H1.5|l and (|3.12() it 
follows that 



d 
dxi 



-'^Casimir (a^l = 0, 2:2 = 0) 



1 



-yL,, 



(5.8) 



where from we conclude, that at T = Tc the Casimir 
force is an increasing function of T for a > 1 (see Fig. 
O, and a decreasing function of T when ct < 1 (see Fig. 
ISJ. Therefore, at the critical point the monotonicity of 
the force changes as a function of cr at cr = 1 where we 
have an inflexion point. 



VI. DISCUSSION 

In the current article we consider the behavior of the 
excess finite-size free energy and the Casimir (solvation) 
force in a classical system with leading long range interac- 
tions in the limit n -^ oo of the 0{n) models (i.e. within 
the spherical model). The dimensionality d and the pa- 
rameter controlling the range of the interaction a are 
chosen so, that the hyperscaling is vahd, i.e. a < d < 2a 
is supposed. In this regime the critical exponents depend 
on a. We demonstrate that, despite of the choice of cr, 
the excess free energy scaling function Xf (see FIG. 01 
and FIG. El is a monotonic function of the temperature 
T and the magnetic field H, with Xf being always a neg- 
ative function. Surprisingly, to a given extend, the above 
properties do not hold in such a general fashion for the 
Gasimir (solvation) force (see FIG. Eland FIG.[7|l. This is 
in line with the results of Section IVJ where we show ana- 
lytically that the force is attractive for any T and cr > 1, 
as well as for any T > Tc ii a < 1. The monotonicity of 
the force turns out to depend on a. For example, if cr > 1 



at T = Tc and H = the force is an increasing function 
of T and L~^, while for ct < 1 it is a decreasing function 
of both T and L^^ at this point (see Eq. H5.8fl and FIG. 
j^J. In addition, one derives that for T = Tc the mini- 
mum of the force is not a.t H = (see FIG. [TJ. Indeed, 
at T ~ Tc the minimum has been found to be at some 
finite value of the scaling field variable X2 ^ HL^/'^ . For 
a = 2, 1.5, 1 and 0.75 the minimum at T — Tc is found 
to be at X2 ^ 0.084, 0.145, 0.263 and 0.416, respectively. 
Such an occurrence of a force minimum for a nonzero 
bulk field has also been reported for the case of (-I-, +) 
boundary conditions |53l l54l|. Here, in this Section, we 
provide more details for the universal finite-size scaling 
function of the Casimir force Xcasimir(a;i, 2:2) presenting 
the numerical results for it as a function of both xi and 
X2 > in FIG. (HI There the effects due to both the tem- 
perature and the magnetic field are demonstrated (we 
recall that xi - (T - Tc)L^^'', X2 - HL^/"). We ob- 
serve that for T < Tc and iJ ^ a cavity shows up in 
the vicinity of the critical temperature that disappears 
for temperatures far away from the critical point. More 
precisely, one observes that there exists a finite value a;* 
of xi, such that for any x* > xi > there is a local min- 
imum of the force at some finite a;2.min, i-C at H ^ 0. 
For xi > xl there is no such minimum at nonzero H. In 
FIG.JHlthe last is shown for the cases a = 0.75, 1, 1.5 and 
cr = 2 (the short-range case). Note, that for a = 0.75 
one needs to go deeply in the under critical region to find 
out where exactly the cavity vanishes. In the short-range 
case cr = 2 we established that xl ~ 0.28. 
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APPENDIX A: SOME PROPERTIES OF THE 
MITTAG-LEFFLER TYPE FUNCTIONS 

The Mittag-leffier type functions are defined by the 
power series |42 |: 



EcA^) = J2 



k=0 



r{ak + (3)' 



a, (3 > 0. 



(Al) 



They are entire functions of finite order of growth. The 
functions are named after Mittag-Leffler who first consid- 
ered the particular case /3 = 1. These function are very 
popular in the field of fractional calculus (for a recent 
review see Ref. |42|'). 

One of the most useful property of these functions is 
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the identity [42 



dxe-'^x'^-^Ec^^p (-x"z) , 



(A2) 



which is obtained by means of term-by-term integration 
of the series (jAlll . The integral in Eq. (|A2p converges 
in the complex plane to the left of the line Re z^/" = 1, 
|argz| < ^a-K. The identity (|A2() lies in the basis of 
the mathematical investigation of finite-size scaling in the 
spherical model with algebraically decaying long-range 
interaction (sec Rcf. [3| and references therein). 

In some particular cases the functions Ea^p^z) reduce 
to known functions. For example, in the case correspond- 
ing to the short range case we have 



Ei,iiz) = exp(z). 



(A3) 



Setting z ~ y ", y > 0, and x ~ ty, we obtain the 
Laplace transform 
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dte-yH^-^Eo,,p {~e) 



1 + 2/° JO 
from which we derive the usefull identity 
1 



(A4) 



1 + z" 



da;exp(-a;z)x" Ea^ai.-x"), (A5) 



by setting (3 = a. 

The asymptotic behavior for large arguments of the 
Mittag-Leffler functions is given by the Lemma |55l| : 

Let < a < 2, /3 be an arbitrary complex number, and 
7 be a real number obeying the condition 



1 



m: < 7 < minJTr, air}. 



Then for any integer p > 1 the following asymptotic ex- 
pressions hold when \z\ —^ oo: 



• At I argzj < 7, 



At 7 < I argz| < tt. 
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APPENDIX B: ASYMPTOTICS OF THE 

FUNCTION /Cd,4j/) 

Here we will evaluate the asymptotic behaviors of the 
auxiliary function lCd,(j{y) used in the expression of the 
free energy and the quantities descending from it. It is 
defined by: 



^d,a{y) 



2(47r)'i/2 7o 
■xE^,i {-x^y 



dxx 



A 



1 
(Bla) 



where 



Aiu)= Y. ^""''• 



(Bib) 



Using the identity 

Ea,l{^^) = 1 ^ zEa,a+li-z), 



(B2) 



it is possible to write down Eq. HBla|l in a more conve- 
nient form, which will allow us to extract the asymptotics 
of the function under investigation. After some algebra 
one obtains 



iCd.Av) = '^^"'/'r (- j ad) - ^XdAv). (B3a) 



where we have introduced the auxiliary function 



idAv) = 



iAirr/- JO 



dxx 2 ' 



-^'4^'-^ 
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(B3b) 



Now, setting x = z(27r) ^ and with the help of the 
identity 



Aiu)^.r-Ar- 

V u \ u 



(B4) 



we rewrite equation (jB3bp (after some algebra) in the 
form 



Xd.aiy) = y- 






dxx 2 2 2 



Jd-l)/2 j-oo 



{2ttY 
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The integral in the second term of the right-hand side 
of (|1J5() can be evaluated exactly with the help of the 
identities 



1 



rH 



and 



dtf-^e-""' 



, m/i' 



Jo ^0^ sm(7r^/i/) 

to yield the result 

2(d-l)-ii?d_i,,2/(''-i)/^ 



(B6) 



(B7) 



(B8) 



For the evaluation of the first integral in the right hand 
side of (|1J5|1 . we note that the two terms in the square 
brackets in ljB5(l cannot be integrated separately, since 
they diverge. Nevertheless, it is possible to outwit this 
divergence, by transforming further IBSp by adding and 
subtracting from the function Ea,a+i{z) its asymptotic 
behavior at small arguments, leading, after some algebra, 
to 

2fw^i^-^^"^^-^''^+^^-(^)- (^^^^ 

Here we introduced the notations 



Cd,<T 


— / dxa;^~2~2 
Jo 


oo 1 — ' 




d-1 < a, 






(B9b) 


and 








^ ' 1=1 •'0 


-i^-xl^ 




X 


^f .f +1 [~yj^^ 


xt \ 1 " 




(B9c) 

1 




inr) r[f + i]J 



Collecting the above results, we obtain 



/Crf,.(y)=a^-'^/2r(^)C(rf) 



(d-l)/a _ 



-a{d-l)-^Dd^i^,y 



-ad-^Dd.^y'"/" - ^TZdAy)- 



T,{d-i)/2 n 



V- 



d,a 



(27r)'^ r[cr/2] 

(BlO) 



The constant Cd,cr introduced in Eg. IBQbll is the so 
called Madelung constant (see e.g. |5a.l57l|') 



a 



d,a 



= lim I 2 V 



1=1 



T[{a-d+l)/2,5P] 

l{<j-d+l)/2 



^' iia-d+i)/2 \,d-l<a, 

(Bll) 
where r[a,x] is the incomplete gamma function. It has 
been shown that this constant has a remarkable property 
of symmetry J57j . which relates its values in the case 
d — 1 < cr to those in the case d — 1 > a. On the other 
hand, it has been shown that Cd.a can be expressed in 
terms of the analytic continuation, over d—1 < cr, of (for 
details see |^) 



Cd.a = 27rs 



7—dj 



d — (7 



Cid-a), d-1 >a. (B12) 



Eq. (jBlOp is the general form of the functions JCd,a{y)- 
According to Eqs. (|B11|) and (|B12|) it can be used to 
investigate the critical behavior of the system for any 
dimension less than d. 

For small y the asymptotic behavior of the function 
^d a{y) is easily deduced from equation ljB10|l . It is given 
by' 



- _^r (f ) ad) - ^%^y(''-^)/^ 

^dAv) « I ^r (f ) ad) - 2y [(4^)-/2^r[a/2]] -' (1 - iny), 

, ^;372i (2 j ^w ~ y (277)- ri572T + ^ d-1 y 



< d- 1 < a, 



= d-l, 



0< o- <d-l. 



(B13) 



For large y the asymptotic behavior of the function 
ICd,a-{y) is obtained by substituting the large x behavior 
of the functions Ea.pix) (given in Eq. ^A7\i ) in the defi- 
nition lBla|l . After some calculations one ends up with 



where 



2i+- r[(d + a)/2] 



(B14b) 



K-d.a{y) — ad,ay , 



(B14a) 
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